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1 Introduction 

The spin problem of the nucleoli has been investigated extensively since the discov- 
ery of the EMC spin effect, which indicates that a small fraction of the nucleon spin is 
from the quark helicity [1]. This surprising result leads to extensive study on the gluon- 
and sea quark-polarization. The current parameterizations suggest the nonvanished po- 
larization in gluon. It implies that the spin-correlation exists among polarized partons. 
Determination of polarized parton distribution functions is crucial for understanding the 
spin structure of the nucleon. The knowledge about polarized parton distributions in 
the nucleon is mainly extracted from global analysis of polarized experimental data at 
different values of Q^ based on the linear (spin-dependent) Altarelli-Parisi (AP) equation 
[2]. However, the contributions of spin-correlation among initial polarized partons are 
excluded in the linear evolution equation. Therefore, it is necessary to modify the AP 
equation to include such correlation effect. 

In the unpolarized case, the parton recombination as a leading correlation effect of 
partons in QCD evolution equation has been studied at leading logarithmic (Q^) approxi- 
mation {LLA[Q'^y) and in a whole x region by [3-4]. We call it the recombination evolution 
equation. The parton recombination effect is generally considered as a phenomenon at the 
small X region, where the rapid increase of parton densities can be obviously suppressed, 
i.e., the saturation is predicted [5]. The behavior of polarized parton distributions at 
small X is an important issue because the experimental determination of the first moment 
of polarized gluon distribution depending on it is rather strongly [1,6]. Moreover, the 
parton recombination can also influence the polarized parton distributions at a larger x 
region. For example, the helicity of a fast gluon may be changed after combining with 
a slower gluon, although its momentum is almost unchanged. Therefore, in this paper 
we generalize the recombination equation to the spin-dependent case. This equation can 



work in a whole x region. Thus, the recombination equation can be used to study the 
spin effect. In particular, we show that the polarized recombination effect disappears at 
the small x limit. We also find that different from the linear evolution equation, the cor- 
relation form of polarized partons in the spin-dependent recombination is relevant to the 
structure of the spin-independent evolution kernel. Therefore, this research is also useful 
for understanding the shadowing dynamics in the unpolarized parton distributions. 

The paper is organized as follows. In Section 2 we construct the spin-dependent 
recombination equation using the same way as refs. [3,4]. In Section 3 the polarized 
recombination functions are evolute. Then the discussion and concluding remarks are 
given in Section 4. 



2 Recombination evolution equation 

The elementary process is one-parton splitting to two-partons in the AP equation. At 
very large number densities of partons, for example in the small x region, the wave func- 
tions of partons can overlap. In this case the contributions of two-partons-to-two-partons 
subprocess (i.e., the parton recombination) should be considered in the QCD evolution 
equations. A complete derivation for such equations need to sum the contributions from 
real and interference Feynman diagrams and corresponding virtual diagrams. Gribov, 
Levin and Ryskin in [7] use the AGK cutting rule [9] to count the contributions of in- 
teference diagrams. On the other hand, the real diagrams are calculated at the double 
leading logarithmic approximation [DLLA) in a covariant way (the cut vertex technique) 
by Mueller and Qiu in [8]. This nonlinear evolution is called the GLR-MQ equation. 

Unfortunately, the application of the AGK cutting rule in the GLR-MQ equation was 
questioned since it breaks the evolution kernels [3]. To avoid this disadvantage, the simple 
relations among the relevant high twist amplitudes are derived by using the time ordered 
perturbation theory (TOPT) at the leading logarithmic approximation {LLA{Q'^)) in [3]. 
Then the recombination functions are calculated in a whole x region based on the same 
TOPT- framework in [4]. This is the recombination evolution equation. 

The helicity amplitudes in the recombination equation satisfy the same TOPT-relations 
as in the unpolarized case. Thus, using the arguments in [3] we can directly write down 
the spin-dependent evolution equations which are 

dxG±{x,Q'^) 



aiK 



rflnQ2 
/ dxidx2[Fig.la + Fig.lb + Fig.lc] 

Jx-\ +XO>X 



CJ Jxi+X2>X 

-2 dxidx2[Fig.la + Fig. lb + Fig.lc], (la) 

J(x-i+x-2)/2>x 



Fig.la = fG+G+ (a:i, 0:2, Q'^)Rg+g+^g± {xi, X2, x) 

+ fG+G- (a:i, Xa, Q^)RG+G-^G±ixi,X2, x) 
+ fG-G+{xi,X2,Q'^)RG-G+~.G±{xi,X2,x) 
+ fG-G- (a^l, a;2, Q'^)RG-G-~.G±iXi,X2, x), (16) 

Fig.lb = /g^^^(a;i, xa, Q^)Rg^^^^^G^ (xi, Xs, x) 



+fqi^-f_{Xl,X2, Q )Rqi^-^_^Q^{Xi,X2, 



~'~-'9L9+V-^l'"^2 5 Q )Rq'_q'j^^G±{Xl,X2, 



X] 



X] 



+fql_q'_ (a^i, a;2, Q^)Rql^q^_^G±i^u^2, x) , (Ic) 

Fig.lc = fc+q^ (a;i, 0:2, Q^)i?G+g;^G± (a^i> 2:2, x) 
+fG+qt (2^15 2:2, Q )-RG+gL^G±(^i5 ^2, a:) 
+/G_g;(a^i, a;2, Q^)i?G_g;^G±(a^i, 2:2, a:) 
+/G^gL (^i> a;2, Q^)^G-<7L-G±(a;i, a;2, a;), (Irf) 

dxq!^{x,Q'^) 
dlnQ2 

, rfxirfx2[Fi^.lrf + Fi^.l(e- l)+Fi^.l(e-2)+Fi^.l/] 

-2 / rfxirfx2[Fi^.lrf + Fig.l{e - 1) + Fi^.l(e - 2) + Fig.lf], (le) 

J(a:i+a:2)/2>a: 



a?ir 



Fi^.lrf = /G+G+(a;i, a;2, Q^)RG+G+^q!^ixi,X2, x) 

+ fG+G-{XuX2,Q'^)RG+G.^q^JXi,X2,x) 



+ fG^G+{xi,X2,Q^)RG_G+^q^^ixi,X2,x) 

+fG^G- {Xl, X2, Q'^)RG^G-^q^ {Xi, X2, x) , (1/) 



Fig.l{e- 1) = fgJ■^-J■^{xl,X2,Q^)RgJ^g,;^^g^Jxl,X2,x) 



+ fgl7ii {XUX2, Q'^)RqJ^q3_^q^^ {xi,X2, 



X] 



X] 



+fqiqi{.xi,X2,Q )R^_-i_^^^{xi,X2,x), {Ig) 



Fig.l{e-2) = f^^^k^{xi,X2,Q )R^j^^k^_^^^{xi,X2,x) 



+ fq^qt (^1' ^2, Q )Rqj^qk_^q^{Xl,X2, x) 



+ fgJ^gi{Xl,X2,Q^)RgJ_gk^_g^jXl,X2,X) 



+ fg^_gk_iXl,X2,Q )Rg,_gk_^g^{Xl,X2,X), {ik) 



Fig.lf = fg^^G+ixl, X2, Q^)Rg^^G+-.q'^ixl, X2, x) 

+ /i3!^G_(^l5^2 5 Q jRqi^G -^q^^\Xli X2-, X) 



'^tcf_G+\XliX2iQ )Rqi_G+^q\\XliX2 



X] 



+ fq'_G- (Xl, X2, Q^)Rg,_G^^q'^{Xl,X2, x) . (U) 

The explanations about this equation are summarized as follows, (i) We neglect the 
linear terms of the AP equation in Eq. 1. (ii) The first integral in Eqs. la and le is from 
two-parton-to-two-parton amplitudes and it leads to the positive (antishadowing) effect, 
while the second integral is the contributions of the interference amplitudes between the 
one-parton-to-two-parton and the three-parton-to-two-parton processes, and it is the neg- 
ative (shadowing) effect. The contributions from the virtual diagrams are cancelled each 



other in the nucleoli, (iii) Ra^b„^Ci{xi, X2, x) are the recombination functions defined as 
Fig.l, where m and n are the hehcities of parton A and B, and we shall detail them in 
next Section, (iv) The nonlinear coefficient K depends on the definition of the two parton 
correlation fAmB„{xi,X2, Q^) and the geometric distributions of partons inside the target. 
For simplicity, we take K as a free parameter. 

The parton density is a concept defined at twist-2. The parton correlation function 
is the generalization of the parton density beyond the leading twist and it has not yet 
been determined either theoretically or experimentally. Similar to [4] , we use a toy model 
(model I) to assume that 

{Xi + X2)fA^B„{Xi,X2,Q'^) = XXlX2Am{Xi)Bn{x2)9{l - Xi - X2), (2) 

where the correlator x is regarded as a constant and it can be incorporated into the free 
parameter K. 

Eq.(l) can immediately be simplified using parity conservation in QCD, which leads 
to 

Ra+B+^C'+ = Ra-B-^C-1 Ra+B-^c+ = Ra-B+^C^) 

Ra+B-^C- = RA-B+^C+1 Ra+b+^c- = Ra-B--*c+- (3) 

We define 

A = A+ + A_, AA = A+- A_, (4) 

and A is the unpolarized parton (quark or gluon) density. In the meantime, the following 
definitions are convenient in the summation of Eq. (1). 



Rab^c — 7:[Ra+b+^c+ + Ra+b+^C- + Ra+b^~^c+ + Ra+b^^C-], (5a) 



^Rab->c — ■7^[Ra+b+^c+ — Ra+b+^C- — Ra+b^^c+ + Ra+b^^C-]-, (56) 



^Rab^c = -^[Ra+b+^c+ — Ra+b+^c- + Ra+b-^c+ — Ra+B-^C-]- (5c) 



The result equation is 



^^^(a;, Q^) 
rflnQ2 

-|-^ / dxidx2 ^-^[Eq.6b] 

-2 / dxidx2 [Eq.6b], (6a) 

J{xi+X2)/2>x Xi + X2 



Eq.Qb = G{xi)G{x2)9{l - xi - X2)RGG^Gixi,X2, x) 
+ m g*(xi)g*(x2)6'(l - Xi - a;2)i?,,g.^G(xi, X2, x) 

+ Y^ G{Xi)q\x2)9{l -Xi- X2)RGq^^G{Xl,X2, x), (66) 

i 

dxAG{x,Q^) 



d\nQ^ 

a^K f XX1X2 



/ dxidx2 [Eq.lh] 

JxT+Xn>X Xl + Xo 



Q Jxi+X2>x Xi + X2 

—2 / dxidx2 [Eq.7b], (7a) 

J{xi+X2)/2>x Xi + X2 



Eg.76 = G{xi)AG{x2)9{l - xi - X2)Ai2^(^^(^(xi, X2, x) 
+AG(xi)G'(x2)^(l - xi - X2)ARgg^g{xi^ 2:2, x) 



+ Yl q'{xi)Aq\x2)6{l - xi- X2)AR^^,-,_^^{xi,X2, x) 

i 

+ ^ Aq\xi)q\x2)9{l -xi- X2)ARI!-,_^q{xi,X2, x) 

i 

+ J2 G{Xi)Aq\x2)e{l -Xi- X2)AR^Qg,^(j{xi,X2, x) 

i 

+ J2 AG{xi)q\x2)e{l -xi- X2)AR^Jq,^c{xi,X2, x), (76) 

i 

dxq^{x,Q'^) 



d\nQ^ 

r XX A Xo 

-2 / dxidx2 -^[Eq.m, (8a) 

J{xi+X2)/2>x Xi + X2 

Eq.Sb = G{xi)G{x2)9{l — Xi — X2)RGG->q'{^ly^2,x) 
+q\Xi)q''{x2)0{l - Xi - X2)Rqjqk^qi{xi,X2,x) 
+q'{Xi)G{x2)6{l - Xi - X2)RqrG^qi{Xi,X2,x), (86) 

dxAq\x,Q'^) 



d\nQ^ 



/ dxidx2 [-Eg. 96] 

Jxt+x-)>x X^ -\- Xo 



Q Jxi+X2>x Xi + X2 

-2 / dxidx2 ^^^^'^ \Eq.9b], (9a) 

J{xi+X2)/2>x Xi + X2 



Eq.9b = G{xi)AG{x2)9{l — xi — X2)ARqq_^^^{xi, X2-, x) 

+AG{xi)G{x2)9{l — xi — X2)ARQQ^gi{xi, X2, x) 

+q^{xi)Aq^{x2)9{l -Xi -X2)ARJ^j-j^^,{xi,X2,x) 

+Aq\xi)q\x2)9{l - xi - X2)AR^J-j _^^,{xi, X2, x) 

+q^ {xi)Aq^{x2)9{l - Xi - X2)ARlj^k^q^{xl,X2,x) 



+Aq^{xi)q''{x2)9{l - Xi - X2)ARgjgk^g^{xl,X2,x) 
+q''{xi)AG{x2)9{l - xi - X2)ARgiQ^gi{xi,X2,x) 
+Aq'{xi)G{x2)e{l - xi - X2)ARl{c^g,{xi,X2,x). (96) 



3 Polarized recombination functions 

The parton recombination function for A{pi) + B{p2) -^ CiPs) + D{p4) in Eq. (1) is 
defined as 

2 D J ulj_ 1 3^33^4 I,, 12 7 dlj_ / -I r\\ 

where 1^ = P^ +Py- The matrix Mp^pj^pgp^ describe the subprocess (for example, Fig. 1), 
which are factorized from the deep inelastic scattering amplitudes due to the application 
of TOPT in [3]. Note that the parton indicated by 'x' is on-mass-shell, (although is off- 
energy-shell) in TOPT, therefore, it has determinate helicity. Concretely, the momenta 
of the initial and final partons are parameterized as 



Pl = [Xip, 0,0,Xip], P2 = [X2P,0,0,X2P], 

p2 _|_ p2 p"^ -\- p^ 
VZ = Np+ ^ -,Px,Py,X3p], ;P4 = [X4P+ -^ ^, -p^, -py, X4P], 

P^ + P^ 
h = [{X4 - X2)p + -^ -, -Px, -Py, (Xi - X2)p], 

P^ + P^ 
Ir = [{Xa - X2)p + -^ -, -Px, -Py, {Xa - X2)p]. (H) 

We take the physical axial gauge and the light-like vector n fixes the gauge as n-A = 0, 
A being the gluon field. Now a recombination function for G±G± -^ ql_ in the t-channel 
is computed in following form. 



1 XX 1 ih 'Y ^a */3 fi *v 

where the polarization vector of the collinear initial gluon is 

e± = [0,l,±2,0]/v^. (13) 



pt — / \ XsX4 ^ n n n ^t^ ,*y ^a *li ^k *c, 

^G±G±~*G± — \q/ J~ — -T^'-^M«;x'-"K</''-""w'-"/3AfTe±e± e±fc± fc±fc± 

o color \-l-\ ~r J^2j 

72/2 r 177*12^' v-'-'iJ 

where r, s are the space indices of p, a of ^4; Car^p and Cpxa are triple gluon vertex and 
the polarization vector of final gluon, which has transverse momentum takes 



1 ±i p^±ipy. 

The gluon polarization tensor on the Feynman propagator is summing over helicity 
state since it is off-shell, i.e.. 



d,M) = 9,. - ^J^^^^j^. (16) 

Similarly, we write other recombination functions in the t-chanel 



gi4-9± ^9 color {xi+X2y ^ ^^ 2 ' ^^'"^ 2 '^ 

xTr[j647« hlpi.-^^)]— — ^272 • (17) 



Dt /1\ ^32^4 ^ r . ,1±75. , /I ±75 NT 

^ri n M*i3d^n[i'L)dux[lR),^rs PiPii /ion 

xGa,,pO^Aae±e± ^272 \" ~ ■^^12 J' ^^^) 

'■l'-r IP4I 



^^-i-«- = ^llolorW^^^''^^'^' A7.(^^)] 

v/^ /^ ^a ^*P^P ^*a df^n{l-L)d,yx{lR) , , 

XL/a^pO^Ao-e±e± ^±^± 7272 ' y^^) 

'■L''R 



where the polarization vector for the on-shell final state gluon (^4), which has transverse 
components is 



1 ±1 P^±iPy 



x|fr-f|]. (21) 

The massless partons with the parallel momenta can go on-mass-shell simultaneously 
in the collinear case and the collinear singularity may arise in the s-channel in calculations 
of the recombination function. In this case, the collinear time ordered perturbation theory, 
which is developed in [10] provides a safe way to evolute such channel. For example. 



3 X3X4 1 ± 75 ^ ^ ^ n'^nP n^n" 



G±G^^q^^ h6^o«or(a;i+ 2:2)3 ^^'"'^ 2 ' "^^'"^h ■ nf {Ir- nf 

xC^^^C'^^^e^ei^e^ef, (23) 



^'^^''^-''^ ~^W color {x,+X,f{h-ny^ ^ W^' 

X tle*ltltie''±e*^ [6^' - p^] (24) 

IP4I 






Ds /2v 3^33^4 ^ ,, 1±75 T^n 1±75 7-n 



xe>'±t*^[6'' - -^l, (26) 



Ds /2, X:iXi 7-n l±75 7-n 



xe^e^'^e^ef, (27) 



and 






V±54-G± '3',oior{x^+X2y "^""^ 2 ''^"'"^ 2 '\lL-nY{lR-n) 



2 



|p4r 



The complete polarized recombination functions, which summing the contributions of 
the t-, M-, s- and their interference channels present in Appendix. At the small x limit, one 
can keep only the InQ^ ln(l/a;) factor (i.e., the DLLA ) and the gluonic initial partons, 
the nonvanished polarized recombination functions can be simplified as 



i?G,G,^G, = ^i(| + | + l), (29) 



and 



Rg+G-^G+ — Rg+G--'G- — 7 -, 72 



9 5 Xi^ + 6 a;2a;i^ + 3 a;2^a;i^ — 4 a;2^a;i^ + 3 a;2^a;i^ + 6 3:2^X1 + 5 X2^ 

4 X {Xi + X2) X'2^X\^ 

(31) 



Further assuming x\ = X2, we have 



271 
Rg+g+-.g+ = Rg+g+^G- = -J-, (32) 



271 
Rg+G-^g+ = Rq+G-^G- = -^ — • (33) 

In this case, Eqs. (5b), (5c) and (7) are zero. It implies that the polarized recombination 
effect disappears. 

The corresponding unpolarized gluon recombination functions at the DLLA are 



9 6 xi" + 8 X2a;i^ + 5 a;2^a;i'^ — 2 a;2^a;i^ + 5 a;2'^a;i^ + 8 a;2^a;i + 6 X2 

4 (xi + X2) X'2^X\^X 



Rgg^g = ; : ^2 — n — ^ (34) 



for Xi 7^ X2, and 



81 1 
Rgg^g — ~r~-, (35) 

4 X 



for xi = X2- 



4 Discussions and summary 

Comparing with the hnear spin-dependent AP equation, the nonhnear recombination 
equation has more comphcated structure, in particulary this equation contains yet un- 
known twist-4 matrix, which is defined as the correlation function fA„^B„ixi,X2,Q^) in 
Eq. (1). The concrete form about the correlation function is an unresolved problem. 
Although we used a toy model (Eq. (2)), we still have other choices. For example. 

Model II. 

fA±,B^{Xl,X2,Q'^) :^ fA±,B±iXi,X2,Q'^), (36) 

i.e., the recombination is dominated by two polarized partons with opposite helicity. In 
this case, we have 

D^Rab-^CD = ^^^AB-^CD^ (37) 

and 

Rab-^cd -^ -[Ra+b^~*c+d + Ra^b+^C-d]- (38) 

Model III. 

fA±,B±iXi,X2,Q'^) > fA±,B^iXi,X2,Q'^), (39) 

i.e., the recombination partons have same helicity. In this case, we have 

^R\b-.cd = ^Rab^cd^ (40) 

and 



1 
Rab^cd —^ -^[Ra+b+^c+d + Ra+b+^C-d]- (41) 



One can find the singularities in the recombination equation (6)- (7), which arise from 
soft initial partons. Since soft partons between the perturbative and nonperturbative parts 
of a DIS amplitude should be absorbed into the nonperturbative correlation function in the 
coUinear factorization scheme, we further modify the above mentioned correlation models 
to avoid the contributions from soft initial partons. A such correlator between two initial 
partons with different values of x is used in [11]. In concrete, we assume that a parton 
with rapidity y = ln(l/a:i) combines with an other parton with rapidity y + r] = ln(l/a;2), 
where X2 = ^Xi and ^ = e~^, one can generalized the correlation function Eq. (2) to 

{Xi+X2)fA,^B„{Xi,X2,Q^) = x{\v\)xiX2Am{Xi)Bn{x2)9{l - Xi - X2), (42) 

where x(l^l) is a normalized correlator of fusing partons with the rapidity difference 77. 
For example, we can take 

X{\v\) = -^e~^^\ (43) 

If a = 5, it implies that \x2 — xi\max = 0.0067. 

In summary, the corrections of parton recombination to the polarized AP evolution 
equation are derived. The relation between spin-dependent and spin-independent recom- 
bination functions and corresponding evolution equations are discussed. 

Acknowledgments: We would like to thank Y.D. Li and J.F. Yang for discussions. 
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Appendix: 

For completeness we list here the results for the polarized recombination functions as 
calculated using equations in Section. 3. 



G^G^ ^G+ = ^_^^l±^^^_(xi^ - 2 x,'x + x,^x' + xs" - 2 xs^x + X2 V 

4:XX2'^ {Xi + X2) X\^ 



,222 2, 2\ /Ai\ 

+Xi X2 — X\ X2X — X\X2 X + X1X2X ) (A-lj 

G+G+ -^ G_ = 3 (6xi^X2X + 6xi^X2^x — 3xi^X2X^ — 7xi^X2X^ 

4xX2^ {X\ + X2) Xi^ 
111 222i/^4 1/^32 o3 2 ^72 3io 4 

+11 Xi X2 X +0X2 XiX + o X2 Xi X — 3 X2 XiX — 7 X2 XiX + 2 X1X2X 

+Xi^ + X2® + 2 Xi^X2 + 2 Xi'^X2^ - Xi V - 2 Xi V + 2 Xi V + 2 X2^Xi 
+2 X2^Xi2 - X2 V - 2 X2 V + 2 X2 V + 2 Xi^X2^) (A.2) 

G+G_ ^ G+ = (a^i + 3^2 - :^) — ^g^ a,^2^4^^4 ^ 282 Xi3x2^x2 + 110 xi Vx2 + 84 xi^X2^x 

8 XX2^ (xi + X2) Xi^ 

— 167xi X2 X — 106x2 xix + 76xi X X2 + 298 xi x X2 — 167x2 xi x 
-310 xi Vx2^ - 106 xi^X2X - 222 xi^x^X2^ + 34 xi^x*^X2 - 75 xi^X2^x 
+218X1^X^X2^ + 110X2^X1X2 - 310xiVx2^ + 81xiVx22 + 84xi^X2^X 
-75xi^X2^x + 282xiVx2^ + 218x1^x2^ - 90xiVx2 - 90x2^Xix3 
+34x2Vxi - 222x1^X2^ + 10xi^° + 10x2^° + 128xi^X2^ + 68xiW 
+24x1^X2^ + 68x1^x2^ + 52x1^X2 + 114x1^X2^ + 25xiV - 24xi^x 
-18xiV + 7xiV + 128x2^x1^ + 114x2^X1^ + 25x2^x2 - I8X2V + 7x2^x^ 
+52x2V-24x2^x) (A.3) 

G+G^ ^G^ = l i^i+^2-x) (246xiVx2^ - 183xi^xx2' - 272xiVx2' 

8 XX2^ (Xi + X2) Xi^ 

+122x2Va;2 + 36xiVx2 - 69xi^X2^x + 107xi^X2^x + 97xi^xx2^ 
-87xi''xx2^ - 103x1^x^X2 + 68x1^x^^x2^ - 105x2^xix^ - 225x1^x^X2^ 
-292 xi^x^X2^ + 240 Xi2x2^x2 + 85 xi\^X2^ + 40 X2^x^xi + 248 Xi^x2x22 



-247x1^X2^ + 77a;iVx2^ + 126xiVx2 + 270x1^x^X2^ - lllxi^xx2 



3^ 5^2 



8^ 2 



7^ 3 



173x2'xi X - 109x2*xix + 266xi''x2*x^ + 114xi*X2^ + 128x/x2 



+68x1^x2^ - 24xi^x + 24x1^X2^ + 25xiV + 68x1^x2^ - 18xiV 



7^ 3 



8^ 2 



7^3 



8^2 



+7xiV + 52x1^X2 + 128x2'xi-' + 114x2 ""xi^ - 18x2'x'' + 25x2 ""x 



(A.4) 



+ 7X2 V + 52x2V - 24X2^X + 10Xi^° + 10X2^°) 

qiqi ^G+ = — -^^^^i^^^^^^7^(36xiVx2 + 8x2^-45xiV^x-9xiX2^x + 72xW2 
27x2 {xi + X2) xxi 

— 36x Xi X2 + 8x1 — 55xi XX2 — 87xi XX2 — 100 Xi XX2 — 88x1 xx2 



3^2^ 3 



-22xi°x^X2 - 60x1^x^X2^ - 22xi''x^X2^ + 8x1^x^x2^ + 222xi^X2 



+256 xi V^ + 222 xi V^ + 136 xi^X2^ + 50 xi ^ - 16 x/x + 8 xi V 



qiqi -^ G- 



+50x2^xi + 136x1^X2^ 

16 (xi + X2 — x) 
27x2 {xi + X2) xxi 



(A.5) 



X2^ — 18x1^x^x2 — 60xi^X2^x — 46xiX2^x + 72x^x1^x2^ 



-90x^xiV^ - 54x^X1^X2^ + 8x1^ - 50xi'^xx2 - 120xi^xx2^ - 138xi'^xx2^ 



5^2 



4^2^ 2 



3^2^ 3 



2^2^ 4 



90 xi XX2 + 43 xi x^X2 + 140 Xi x^X2 + 156 xi x^X2 + 52 xi x^X2 



-96 xi^X2^ + 112 xi V^ + 96 xi V^ + 64 xi^X2^ + 32 xi ^ - 8 xi^x + 4 xi V 



3^ 4 



+32x2^x1 + 64xi X2^ + 5x2^x^x1 + 18xix'^X2** - 16x2^x + 8 x2V) (A. 6) 



qiqi. -^ G+ 
qiqL -^ G- 
G+q\ -^ G^ 







2 (xi + X2 — x) 

9x2X (Xi + X2) Xi^ 



(A.7) 
(A.8) 



3^ 3 



(13xi X2^ - 27xi^x + 27xiV - 5xiV + 16x2-^x1 



4_ 2 



4^2 



3^3 



+24 X2 xi + 16 X2 xi — 3 X2 X — 6 X2 X + 5 X2 x — 45 xi xx2 — 30 Xi XX2 



+36 xi Vx2 + 20 xi Vx2^ - 5 xi^xV + 9 xi® + 4 X2® + 6 X2 Vxi + 18 xi ^ 



-35x2 xi X — 26x2 Xix + 4x2 xix 1 



G+g^ -^ G_ 



1 (Xi + X2 - x) 

18 X2X (xi + X2) xi' 



3^2 



(A.9) 



„3^ 2 



(36 xi - 55 X2 V - 71 X2 X + 110 X2 V - 9 x-^xi + 16 X2 



3^2 



+54 xi XX2 — 88 X2 xix — 9 xi X X2 — 27 xi X + 37 X2 xi x + 48 X2 xi 



+48 a;2^xi + 36 xi^xs + 16 0:2^X1^ + 56 xa^xix^) (A.IO) 

G^q[_ -^ G^ = 3 (8 XiX2"^ — 14 Xia;x2^ + 5 0:1X^X2 + 8 X2^ — 40 X2^x + 64 0:^X2^ 

Q^i^a; (xi + X2) 

— 32a;'^a;2 — 18x1^x2 — 18x1^x2^ + 53xi^xx2 + 27xi^x — 36xi^x^ + 9x^xi) 

(A.ll) 

G+gi -^ G- = — '^^i+^^-x) ^^g^^5^ _ 23a;iV + 88x/xx2 - 36x/x2^ + 135xi^xx2^ 
18 X2X (xi + X2) xi^ 

+7xiV - 54xiVx2 - 72x2^X1^ - 20x2''xi^ - 25xi^x^X2 - 13xi^x^X2^ 
+49x2^xi^x — 7x2^x^xi + 32x2^xi — 53x2'^xix + 28x2^xix^ + 16x2^ 
+30x2V-7x2V-39x2^x) (A.12) 

G+G+ ^ q+ = 77-^—7 — — — ^ — (2x1"^ + 3x1^X2 -4xi^x + 2xi\2^-2xi^X2X + 4x\i^ 

6x2^ (Xi +X2) Xi^ 
+3X2^Xi — 2X2^XiX + 2X2*^ — 4X2^X + 4x^X2^ — XiX2X^) (A. 13) 

G+G+ ^ qt = — — — — 1^ — (4x^xi^ + 4x^X2^ - xiX2X^ + 4xi^X2^ + 4x2^Xi - 4x2^XiX 

6x2^ (Xi + X2) Xi^ 

+2x2^ -4x2^x + 2xi^ + 4x1^X2 -4xi^x -4xi^X2x) (A.14) 

G+G_ -^ q^^ = 7 (6 Xi'^xx2^ — 61 Xi^xx2^ + 6 Xi^xx2^ + 48 x^X2^xi 

12 X2^ (Xi + X2) Xi^ 

+66 xi*^x^X2^ + 66 xi^x^X2^ — 61 xx2^xi^ + 48 xi^x^X2 — 49 X2^xix 
-20 xi Vx2^ - 49 xi^xx2 + 4 xi® + 4 X2® + 16 x/x2 + 24 xi^X2^ + 16 xi^X2^ 
+8 xi*^X2^ - 8 xi^x + 8 xi V + 16 X2^xi^ + 24 X2^xi^ + 16 X2'^xi - 8 X2'^x 
+8xW) (A.15) 

r^ r^ % V"^! + X2 — Xj , 242 on 4 3irr224 01 6 

G+G_ ^^ q^ = 7 (,66 X X2 Xi — 30 X\ xx^ + 66 x X2 Xi — 31 Xi XX2 

12 Xi^ (Xi + X2) X2^ 

+48 Xi^x^X2 + 48 x^X2^Xi — 30 Xi^xx2*^ — 43 xxi^X2^ + 8 x^X2^ — 31 X2^XiX 
-20 x^X2^xi^ - 43 xi^xx2^ - 8 X2^x + 4 xi* + 4 X2^ + 16 xi'^X2 + 24 xi^X2^ 
+16 xi^X2^ + 8 xi*^X2'^ - 8 x/x + 8 x^x^ + 16 X2^xi^ + 24 X2^xi^ 
+16x2^x1) (A.16) 



qW+ ^ql = ^ (^1 + ^2 x)x^ ^^ a;i^ + 12 Xi^X2 + 18 ^1^2' + 12 XiXs^ + 3 Xs^ + 24 Xa'x" 

2/ X2 (xi + 3:2) 

+4:XX2Xi'^ + 8xX2'^Xi + 4:XX2'^) (A. 17) 

• • . 64 a;ia;2(a;i + a;2 -x)^x2 / a 1 o\ 

g^gi ^ gL = — ^- -7 (A. 18) 

y [xi + X2) 

qiqt ^qi = (A. 19) 

qiqi -^ gi = (A.20) 

J ^k , ^» _ 8 (a^i + 3^2 - xf {xi^ + X2^) . . 

y X2 [Xi + 3:2) Xi 

qiqX ^ql = (A.22) 

7 A- » 8 (2^1 + a;2 - x) xi , . ^^, 

qiql - ^V = q ^ . ^ ^3 (A.23) 

y a;2 [Xi + a;2) 



^j ^k _^ J _ 8 (a:i + 0:2 - x) X2 

9 Xi(xi+X2) 



g^g_ -> gL = - ^ ^ — -^ (A.24) 



g^G+ ^ g+ = — ^^^^^^ ^> (128x1^X2^ + 32 xi^ - 126x2^xix + IO8X2VX1 + 86xiV 

18 X2^ (Xi + X2) XiX 

-35 xi^x^ + 72 X2^ - 108 X2^x + 104 Xi2x2^ + 144 X2^Xi + 81 X2^x2 - 29 x^X2^ 
+128x1^x2 + 192x1^X2^ - 83xi^x - 29xWa:i - 263xi^xx2 + 182xiVx2 
-191xi^xx2^ - 29xi2xx2^ + 127x1^x2x2^ - 31x1^x^X2) (A.25) 

q!^G+ ^ gi = (A.26) 

g^G. ^ g+ = — '^^i+^^-a:) ^^^^^^ _ ^^g^^s^ ^ 64xi^X2 - 40xi^X2^ - 143xi^X2X 

18 XXi (Xi + X2) X2'^ 

+158xi'^x2 + 145xi^X22x - 71xi^x^ + 148xi^X2x2 - 144xi^X2^ - 53x1^x^X2 
+295xi2x2^x - 93xiVx22 - 72x12x2"^ - I26X1X2V + 142xiX2^x 
-7 xix3x22 - 23 X2 V + 16 X2^x + 7 x^X2^) (A.27) 

q^G. ^ gi = (A.28) 



GG^G= — — ^^] ^ "^^ ^^^ — (900 Xi^x2x2^ - 408 xi^xxs^ - 882 xi Vx2^ + 284 X2^xix2 
16 XX22 (xi + X2) Xi^ 

+ 100 xi^x^X2 - 198 xi^X2^x + 153 x/x2^x + 143 Xi^xx2^ - 216 xi^xxs^ 



6^3, 



3^4^ 3 



6^ ^3 



5^3^ 2 



3^3^ 4 



265 xi'^x'^xs + 228x1^x^X2'' - 267x2''xix'' - 639xiVx2^ - 902xi''x''x2 



2^4^ 4 



5^4, 



6^2^ 2 



2^ 5^3 



606 Xi^X2*'x^ + 232 xi^x*X2* + 104 X2^x*Xi + 614 xi*'x^X2^ - 661 Xi^X2^x' 



+224 xi Vx2^ + 288 xi Vx2 + 836 xi Vx2^ - 251 Xi^xx2 - 398 X2^xi2x 



3^ 5^2 



8^ 2 



7^ 3 



-249 X2"xix + 832 xi^X2^x^ + 292 xi«X2^ + 360 Xi 'x2'' + 260 Xi V2 



8^2 



56 xi^x + 176 Xi^X2^ + 60 Xi V + 260 Xi^X2'' - 48 Xi ' x'' + 20 Xi V 



8^ 2 



IJi 



8^2 



6^4 



+128x1^x2 + 360x2'xi^ + 292x2*xi^ - 48x2'x^ + 60x2*x^ + 2OX2V 



+128 X2V - 56x2^x + 24xi^° + 24x2^°) 



(A.29) 



j-j n ^6 (X1 + X2-X) 43 

gJgJ ^G = — ■ -^(18 Xi x'^ 

^( XX2 (Xi + X2J Xi 



2^ 5„ 



X2 + 8 X2* - 84 Xi^X2''x - 32 XiX2^X + 72 X*^Xi^X2^ 



-18x^X1^2^ - 18x^xiW + 8x1^ - 32xi^xx2 - 84xi^xx2^ - 148xi^xx2^ 

-148 Xi^XX2'^ - 11 Xi^X^X2 + 10 Xi'^X^X2^ + 50 Xi^X^X2^ + 10 Xi^X^X2^ 

+222 xi^X2^ + 256 xi^X2^ + 222 xi^X2^ + 136 Xi2x2^ + 50 xi^X2 - 8 xi^x 



+4 xi^x^ + 50 X2 xi + 136 xi^X2^ — 11 X2^x^xi + 18 xix^X2*^ — 8 X2 



X 



Gq' ^G 



+4X2V) 
1 (xi + X2 — x) 



(A.30) 



3_ 3 



(16 xi X2^ - 46 xi^x + 29 xi V - 1 1 xi V - 8 X2''xi 



I8X2X (Xi + X2) Xi^ 

+76 X2^xi^ + 96 X2^xi - 101 X2^x + 122 X2 V - 53 X2 V + 8 Xi*^xx2 



3^2 



2^2^ 2 



2^3, 



+133 Xi XX2 — 9 Xi X X2 + 26 Xi X X2 — 18 xi x X2 + 36 Xi + 32 X2 



GG^q' 



-42x2Vxi + 72x1^X2 - 32x2''xi^x - 212x2^xix + 174x2''xix^)(A.31) 
1 (X1+X2-X) qi42,2^^4^^2_^^2xi^xx2=^ + 114x2x2W-78xi6xx2 



12 Xi^ (Xi + X2) X2^ 

+78x1^x^X2 + 78x^X2^X1 - 112xi^xx2'^ - 130xxi\2^ + 16x^2*^ 



-78 X2^xix + 32 x^X2^xi^ - 130 Xi^xx2^ - 16 X2^x + 8 xi^ + 8 X2^ + 39 xi^X2 



5^ 3 



+82xi'^X2^ + 105xi°X2-' + 108x1^x2^ - 16xi'x + 16xi'^x^ + 105x2^xr 



+82X2*^X12 + 39X2^X1) 



(A.32) 



q^q^ ^qi = ^^^ + ^^ — ^±^ (3 a;/ + 12 xi^Xs + 18 Xi^a^ + 12 ^1X2^ + 3 X2^ + 48 ^2 V 
27 a;2 (xi + X2) 

+4xx2Xi^ + 8xx2^Xi + 4xx2^) (A. 33) 

J fc . _ 8 (xi + a;2 - x) (xi^ + a;2^) 
9 0:2 {xi + ^2) xi 



g^g'^ -.q^ = JJ IJl^jtlZ^L^iL^:^ (A.34) 



g^G ^ g* = — '^^^^^^ f^-* (32x1^ - 9X1X2V - 23xi^X2^x + 8xiX2^x - 18xixW 

18xXi (Xi + X2) X2^ 

+ 165xi^X2X^ — 203xi'^X2X — 42xi^x^X2 + 133xi^X2^x + 72x2^Xi 

+ 17X1^X^X2^ + 29 x2V -46X2^X- 11X^X2^- 101 Xi^X + 96X1^X2 

+76 xi^X2^ + 122 xi V - 53 xi^x^ - 8 xi^X2^ + 16 Xi2x2^ + 36 X2®) (A.35) 

they consist with the definition Eq. (5a). 
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Figure Captions 

Fig. 1 The diagrams contributing to the polarized recombination functions. Here the 
shaded part imphes all possible QCD-channels and x means the probing place. 
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